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Introduction 


When the source and observation points are close, the potential integrals over 
wire segments involving the wire kernel are split into parts to handle the singular 
behavior of the integrand [1]. The singularity characteristics of the gradient of the 
wire kernel are different than those of the wire kernel, and the axial and radial 
components have different singularities. The characteristics of the gradient of the 
wire kernel are discussed in [2]. 


Integration of the Gradient of the Wire Kernel 


To evaluate the near electric and magnetic fields of a wire, the integration of the 
gradient of the wire kernel needs to be calculated over the source wire. Since the 
vector bases for current have constant direction on linear wire segments, these 
integrals reduce to integrals of the form 
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where 5 (£j,£ 2 ) is a basis function in nonnalized segment variables c, , =l-£i 

for scalar current, VK(p,z- z') is the gradient of the wire kernel, z is the 
projection of the observation point onto the wire, p is the perpendicular distance 
from the wire axis to the observation point, and Az, and Az 2 are defined as in [1]. 


On the wire surface, the gradient of the wire kernel VK(p,z - z') is dominated by 
a near-delta function behavior in a neighborhood ( Az, and A z 2 ) about the 

projection of the observation point onto the segment (^j 0 ,£ 2 j [2], The integrals L 

and h are designed to account for this behavior. Outside this neighborhood, the 
gradient of the wire kernel has a l/ IV behavior that is treated in I\ and I 4 . 

The first integral, 
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has an integrand dominated by a 1 / R 3 behavior. This behavior is canceled by 
letting 
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or, upon integrating, 
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The use of (3) and (4) allows (2) to be rewritten as 
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and Ac = a/AT in this work. The integral I 4 is similarly evaluated. 

The near-delta function behavior (p near but not equal to a) present in the second 
integral h is accounted for by letting 
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which yields, after integrating, 
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and 


z = z + (p-a) tan [(/?-«)«]. 
Hence, h may now be written as 
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The integral / 3 is evaluated using the same approach as h. 
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The sample points qf* are given in (9) and (17), while the bracketed quantities in 

(5) and (13), represent the new weights. These may be used in typical numerical 
quadrature routines. 


Results 

The near field of a 1 -meter wire dipole with a radius of 1 mm is calculated and 
compared with results using EIGER [3]. The dipole is modeled using 20 linear 
segments and lies on the z axis. It is exited at the center with a unit-strength 
voltage source at 300 MHz. The magnetic field is sampled 1 mm from the dipole 
surface and is shown in Figs. 1 and 2. There is good agreement between EIGER 
and the data generated using this fonnulation. 

Summary and Conclusions 

A fonnulation for integrating the gradient of the thin wire kernel is presented. 
This approach employs a new expression for the gradient of the thin wire kernel 
derived from a recent technique for numerically evaluating the exact thin wire 
kernel. This approach should provide essentially arbitrary accuracy and may be 
used with higher-order elements and basis functions using the procedure 
described in [4]. 
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dipole. 



Fig. 2. The phase of they component of the magnetic field 1 mm from the wire dipole. 



